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ABSTRACT 

We have detected a strong periodicity of 1.80±0.05 hours in photometric observa- 
tions of the brown dwarf Kelu-1. The peak-to-peak amphtudc of the variation is 
~1.1% (11.9±0.8mniag) in a 41nm wide filter centred on 857nin and including the 
dust/temperature sensitive TiO & CrH bands. We have identified two plausible causes 
of variability: surface features rotating into- and out-of-view and so modulating the 
light curve at the rotation period; or, elliposidal variability caused by an orbiting com- 
panion. In the first scenario, we combine the observed vsini of Kelu-1 and standard 
model radius to determine that the axis of rotation is inclined at 65±12° to the line 
of sight. 

Key words: Stars: brown dwarfs - Stars: oscillations - Stars: rotation - Stars: 
atmospheres ~ Binaries: close 



1 INTRODUCTION 

The study of rotation and variability in main sequence stars 
has led to a great improvement in our understanding of 
their physics (e.g. Stauffer & Hartmann 1986). Recently, 
several groups have shown that variability can also be de- 
tected in substellar brown dwarfs (Tinney & ToUey 1999; 
Bailer- Jones & Mundt 1999, 2001; Martin, Zapatero Osorio 
& Lehto 2001). 

In this paper we present differential photometry of the 
brown dwarf KeIu-1, in a search for rotational variability. 
Kelu-1 is a field brown dwarf, discovered by Ruiz, Leggett 
& AUard (1997) via its large proper motion. It is classi- 
fied cis an L2 dwarf in the scheme of Kirkpatrick et al. 
(2000), and model fits to its spectrum estimate an effec- 
tive temperature of 1900±100K (Ruiz et al. 1997). The par- 
allax distance of 19.6±0.5 pc gives an absolute magnitude 
of Mj = 11.96±0.09 (Kirkpatrick et al. 2000) or a bolomet- 
ric magnitude of M)-,q[ = 13.9±0.1. The detection of strong 
lithium absorption (Ruiz et al. 1997) in an object of this lu- 
minosity implies a mass below O.Q7M0 (Tinney 1998), mak- 
ing Kelu-1 a bona fide brown dwarf. In common with many 
L- and M-type brown dwarfs (Basri et al. 2000, Tinney & 
Reid 1998), Kelu-1 is a very rapid rotator with a measured 
ti sin i=60±5 kms~^. It is thought this may indicate that the 
magnetic braking mechanisms which operate in more mas- 
sive stars do not operate with the same efficiency in brown 
dwarfs (as first suggested by Tinney & Reid 1998). 

In section 2 we describe the data acquisition and reduc- 
tion. Time series analysis of the resulting differential pho- 



tometry is presented in section 3. In section 4 we discuss 
our observations in terms of surface features on the brown 
dwarf, or the nature of an orbiting companion. 



2 DATA 

2.1 Observations 

Kelu-1 was observed on 2000 March 24-25 using the Taurus- 
2 instrument on the 3.9m Anglo- Australian Telescope with 
a CCD detector denoted MITLL3. MITLL3 is a 2096x4096 
pixel device with 15^m pixels, giving an image scale of 
0.373" per pixel. Observations were performed through an 
intermediate bandpass blocking filter for the Taurus Tun- 
able Filter. This filter (denoted R6) is 41nm wide and is 
centred at 858nm. The bandpass of the filter is shown over- 
laid on a spectrum of Kelu-1 (taken from Martin et al. 1999) 
in Figure ^ Also plotted in Figure ^ is the spectrum of the 
slightly hotter LO dwarf DENIS-P J0909-0658 (Martin et al. 
1999), indicating the temperature sensitivity of the TiO and 
CrH molecular bandheads selected by our filter. 

A specially constructed slot mask was installed in the fo- 
cal plane wheel of Taurus-2 (the same mask used by Tinney 
& ToUey 1999). Multiple exposures of a sky field can then 
be made through this slot, between which charge is shuffled 
and stored on the unilluminated surface of the CCD. When 
used to observe Kelu-1, the entire Taurus-2 instrument was 
rotated at an angle of 140°, so that several suitablephoto- 
metric reference stars could be observed (see Figure b|) . This 
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Table 1. Log of Kelu-1 observations in March 2000. 
Date Start UT End UT # Images Exposure Time (s) 



each individual exposure by fitting a profile to several of 
the brightest (non-saturated) stars in the field. This profile 
was then scaled to all the stars in the field to determine 
their brightness. No faint neighbours were found after the 
PSF had been subtracted from the target stars. Pixels which 
varied by more than 5% from the mean in the flatfield were 
flagged as bad. These bad pixels were interpolated over dur- 
ing the PSF fitting stage, but completely ignored during 
the actual photometry. PSF photometry also has the ad- 
vantage of being relatively immune to cosmic ray events, 
which would affect aperture photometry. In order to verify 
our PSF photometry, we also performed aperture photome- 
try in the same fashion as Martfn et al. (2001). The aperture 
photometry was found to give the same, but noisier, results 
as the PSF photometry. 

To obtain high precision differential photometry of 
Kelu-1, we compared its apparent brightness with that of 
surrounding reference stars, cancelling the effects of appar- 
ent brightness changes due to variable extinction, seeing, 
instrument performance or exposure time. There are four 
available reference stars in our field-of-view (the brightest 
star in the field was often saturated). We define the bright- 
est of these four stars to be our comparison star, which sets 
the magnitude zero-point for each exposure. The three re- 
maining stars can then be used as check stars, allowing us to 
detect intrinsic variability in our comparison star. One of the 
check stars was found to be variable and therefore discarded. 
The remaining checks were averaged to define a mean check 
magnitude. These stars are identified in Figure ^. 

The differential lightcurves obtained are shown in Fig- 
ure 1^. The top panel shows Kelu-1 minus comparison, and 
the lower panel comparison minus check. It is clear from 
these plots that Kelu-1 minus comparison photometry is 
more variable than comparison minus check photometry - 
indicating the variability detected is due to changes in the 
apparent brightness of Kelu-1 itself (an increase in Sra rep- 
resents a dimming of Kelu-1). Variability in the comparison 
star would manifest itself as equal, but opposite, changes in 
each lightcurve. This effect can be seen on the first group of 
datapoints - the Kelu-1 minus comparison data is above the 
mean, and the comparison minus check photometry is below 
the mean, indicating the comparison star has dimmed. We 
neglect this group of (60) datapoints in further analysis. 

As the photometry presented in Figure ^ is differen- 
tial, most sources of systematic error which might produce 
a "spurious" signal in Kelu-1 are cancelled out. The remain- 
ing possible sources of error are differential effects due to 
the differences in colour of Kelu-1 and the reference stars, 
or their different locations on the CCD. 

Differential position effects have been minimized by 
keeping all objects within 30 pixels of a nominal position on 
the detector for all observations. However, differential pho- 
tometry for Kelu-1 versus detector position shows a residual 
correlation between 5m and y-position on the CCD (Fig- 
ure ^). We have fitted this correlation with a straight line 
and subtracted it. The cause of this effect remains unknown. 
The maximum pixel-to-pixel sensitivity variation in the fiat 
field is 2.5% in the region of interest, and does not display 
any gradient. All data discussed hence have had this corre- 
lation removed. No correlation was found with x-position. 

Differential colour effects will be produced when the ef- 
fective wavelength of a star through a filter is different from 
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Figure 1. The spectrum of KeIu-1 (solid Une) compared with 
the spectrum of the LO (slightly hotter) brown dwarf DENIS-P 
J 0909-0658 (dashed line). The transmission profile of the R6 filter 
is overplotted (dot-dash line) , showing how it isolates the sensitive 
TiO and CrH features at ~850nm. 
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observing technique results in a series of 30 exposures, fol- 
lowing which the CCD is read out only once, providing very 
low read-out overheads and precise timing. Table ^ gives a 
log of the observations. 

2.2 Reduction and photometry 

Each CCD frame contains 30 seperate exposures, so the data 
processing is slightly more complicated than for standard 
CCD imaging. Each frame was bias subtracted and then 
sliced into thirty individual exposures. The header file for 
each exposure was modified so as to include the correct tim- 
ing and positional information. A master flatfield exposure 
was then constructed from dome flats observed in the same 
manner as the data. The data were then divided by this flat, 
providing correction to the ~1% level. 

Due to the proximity of a nearby star to our primary 
target, point-spread function (PSF) photometry of this data 
is required. Photometry was performed using DAOPHOT 
within the IRAF environment. A PSF was determined for 

Figure 2. "Slot" mask image of Kelu-1 and the comparison and 
c/iecfc stars used in our differential photometric analysis. The ver- 
tical bars in the image are bad columns on the M1TLL3 detector. 
The orientation of the sky through the mask when Taurus-2 is 
rotated to PA=140° is indicated. 
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Figure 3. The top panel shows differential photometry f or K elu-1 
(in the sense Kelu-1 minus comparison as described in j |3.2| ) over 
two nights. The bottom panel shows comparison minus check to 
the same scale, clearly indicating the variability detected is from 
Kelu-1, not the comparison. The uncertainties plotted are based 
on the photon counting errors as produced by DAOPHOT and 
propagated to the differential results for Kelu-1, comparison and 
check stars. This plot shows raw differential photometry, before 
the correlation of Figure ^ is removed. 
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Figure 4. Differential photometry plotted against Y pixel po- 
sition on the CCD. There is a clear correlation at the level of 
O.OOSmag per pixel, which we remove by subtracting a straight 
line fit. 
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the reference objects. The target will then suffer a different 
amount of atmospheric extinction as the target rises and 
sets. If such an effect is present we would expect that spuri- 
ous periodicities would be produced at aUases of 24 hours (ie 
6, 12 or 24 hours). We would also expect to observe a cor- 
relation between differential photometry and airmass if this 
effect is present. Figure ^ shows a plot of the differential 
photometry for Kelu-1 (upper panel) and comparison stars 
(lower panel) versus airmass. There is no evidence (from 
Rank-Spearman tests) for a correlation in either of these 
plots. In addition, Figure |^ shows that Kelu-1 does not have 
a strong colour gradient over our filter band-pass. 



3 TIME SERIES ANALYSIS 

We have calculated the least-squares weighted power spec- 
trum of the differential lightcurve to look for periodicities. 
The power spectrum is calculated in the manner described 
by Frandsen et al. (1995), which we briefly describe here. 
The data are represented by the function Ai sm{2n fit + 
for a range of frequencies fi, where Ai, (f)i and t are the 



Figure 5. Differential photometry plotted against airmass for 
Kelu-1 minus comparison (top panel) and comparison minus 
check (bottom panel). A Rank-Spearman correlation test shows 
no significant correlation between airmass and differential mag- 
nitude. 
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amplitude, phase and time respectively. The power at each 
frequency is then Pi — — o? -\- (3^ , where a and {3 are 
given by 



a = 
P = 
with 



(SC2 — Cx)/{S2C2 — X^), 
(CS2 — SX) I (S2C2 — X^) 



C2 
S2 



^ ^ WjXj COs{0,itj), 

y^^WjXjSmjnitj), 
^ ^ Wj cos^ {0,itj), 
^ ^ Wj sin"^ (Qitj) , 
X = y ^ tt)j sin(r2itj) cos(n.itj) 

where {tj,Xj,Wj) are the data points & associated weights 
and Q,i — 27r/i. Each data point is assigned a weight of 
Wj — l/e|, where Cj is the error associated with point j. 
This system gives data points with smaller errors higher 
weights. 

The discrete sampling of observations means that the 
observed power spectrum is the convolution of the "true" 
power spectrum, with the "window function" . We have cal- 
culated the window function as prescribed by Roberts et al. 
(1987); 

W^(/) = ^^e-^-'^'^ (1) 
J 

The power spectrum and window function of our data 
are shown in Figure ^ The upper panel shows the power 
spectrum over the complete range of frequencies available 
to us; one day~^ up to the Nyquist frequency of 0.5 min~^ 
(1 < / < 720 day-i; 24 > P > 0.033 hours). There 
is only one main peak in the power spectrum. The lower 
panel shows a close-up of the region of interest (1 < / < 
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30 day~^). The power spectrum of the Kelu-1 differen- 
tial photometry (solid line) shows a strong peak at a fre- 
quency of 13.34±0.38 day~^, corresponding to a period of 
1.80±0. 05 hours. The peaks at higher and lower frequencies 
are produced by the first sidelobes of the window function, 
as shown in the middle panel of Figure ^ The maximum 
power of the main peak is 3.5xl0~^ rnag^, corresponding to a 
peak-peak amplitude in the lightcurve of 11.9 mmag. Monte- 
Carlo simulations of the data give a 1-cr error in the ampli- 
tude of 11.9±0.8 mmag. CLEANed power spectra (Roberts 
et al. 1987) do not reveal any other periodicities hiding in 
the sidelobes of the main frequency. 

Also shown in the lower panel of Figure ^ is the 
power spectrum of the comparison minus check photome- 
try (dashed line). There is no peak at the frequency de- 
tected in the Kelu-1 photometry, and the maximum power 
in this spectrum is an order of magnitude lower. This con- 
firms the periodic variability we detect in the Kelu-1 mi- 
nus comparison photometry is due to intrinsic variability of 
Kelu-1 rather than the comparison star. 

To further investigate the periodicity detected, we have 
treated each night's data separately. The periodograms are 
shown in Figure |^ The same ~1.8 hour periodicity is de- 
tected in both datasets, indicating the source of the vari- 
ability is stable over at least ~48 hours (~25 rotations). 

Close examination of the data reveals a slow linear 
change in Kelu-l's brightness at a level of 0.003 mag/day 
over the duration of the observations (this is not caused 
by the previously discussed correlation between Sm and y- 
position on the CCD (Figure ^ . We have not removed this 
effect for the time-series analysis. Figure ^ shows the Kelu- 
1 photometry folded on a period of 1.80 hours (with the 
slow, linear trend this time removed), and combined into 5 
minute bins. We derive the following ephemeris for the time 
of minimum light, tj^jj^; 

Vin = JD 2451628.037 + 0.075 x e 
where e is the number of rotations since JD=2451628.037. 



4 DISCUSSION 

Kelu-1 clearly displays periodic variability, but what is the 
cause of this variability? We have developed four possible 
explanations which we discuss in this section: 

(i) surface inhomegeneity moderated by meteorology and 
variable dust formation; 

(ii) surface inhomegeneity moderated by magnetic 
starspots; 

(iii) light curve variability due to gravitational (i.e. tidal) 
distortion of Kelu-l's envelope by a close companion; and 

(iv) light curve variability due to an eclipsing binary. 



4.1 Dust cloud meteorology 

At the cool temperature of Kelu-l's photosphere, dust for- 
mation will play an important role. Theoretical spectra in- 
dicate that Kelu-1 is better matched by models with dust 
suspended in the photosphere than models in which con- 
densates are "rained-out" (Baraffe et al. 1998). We can- 
not expect, however, the atmospheres of brown dwarfs to 



Figure 6. The lower panel shows the least-squares weighted 
power spectrum of Kelu-1 minus comparison (solid line) and com- 
parison minus check (dashed line). There is a clear peak at 13.33 
day~^ (period=1.80±0.05 hours) with a peak-to-peak amplitude 
of 11.9 mmag in the Kelu-1 power spectrum. The error bar repre- 
sents the 1-(T width of the central peak. There is no correspond- 
ing peak in the comparison minus check power spectrum, showing 
that the detected period is from Kelu-1. The middle panel shows 
the window function of the dataset, clearly indicating the peak 
is not due to sampling effects. The upper-most panel shows the 
least-squares weighted power spectrum of Kolu-1 over the entire 
range of frequencies we can study, from 1 day~^ to the Nyquist 
frequency of 0.5 min~^. This panel shows there is no significant 
structure in the power spectrum above f~20day~^ 
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be placid, homogenous places — especially if rotating at 
V sm I =60kms-^ They are no doubt dynamic and evolv- 
ing, much like the atmospheres seen on planets in our own 
solar system. This leads us to investigate the possibility that 
we are observing dust cloud meteorology in Kelu-l's atmo- 
sphere. 

Gelino et al. (2001) suggest there is a possible corre- 
lation between J-Ks colour and variability in L dwarfs - 
variable objects tend to have bluer J-Ka colours, which they 
claim is evidence that holes in dust clouds are causing vari- 
ability. Kelu-1 has a J-Ks ~ 1.6 (Kirkpatrick et al. 1999). 
This is not significantly different from the mean colour for an 
L2 dwarf (J-Ks ~ 1.5), so Kelu-l's strong variability sheds 
no light on this suggestion. 

Could the variability be due to a large feature on one 
side of the brown dwarf? The surfaces of all the Solar System 
giant planets are dominated by surface banding, which is 
driven by an interior in which Coriolis forces dominate over 
buoyant convection, to produce long, thin quasi-cylindrical 
cells oriented with their axes parallel to the axis of rota- 
tion (Schubert & Zhang 2000). In brown dwarfs, however, 
- even brown dwarfs with a 1.8h rotation period - buoy- 
ancy forces should dominate to produce chaotic three dimen- 
sional convection. So banding structures, like those seen on 
Jupiter, are probably not to be expected. If cloud formation 
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Figure 7. Power spectrum of data from night 1 (upper panel) & 
2 (lower panel). The same peak at ~13 day~^ is detected in both 
datasets. The upper section in both panels shows the window 
function. 
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Figure 8. Kelu-1 photometry folded onto the main period (1.80 
hours), and combined into 5 minute bins to improve the signal to 
noise ratio. This data has had a long term trend of (5m brightening 
at ~3mmag/day removed. 



is inhomogenous, rotation and convection should combine to 
produce less ordered structures. So convective upwelling, or 
large cloud inhomogeneities may be possible. 



4.2 Magnetic starspots 

Another possibility is that the imhomogeneities we detect 
are caused by magnetically induced starspots, analogous 
to those observed on the Sun and other cool stars. How- 
ever, Gelino et al. (2001) have shown that the atmospheres 
of brown dwarfs are unable to sustain starspots, as low 
Reynolds numbers imply the atmosphere is not bound to 
magnetic field lines. In addition, several groups have shown 
that there no strong evidence for the link between variability 
and activity (as measured by H-a equivelent width), which 
would be expected if magnetic forces drove variability. We 
note however that this conclusion is weakened by the small 
number statistics currently involved. 



4.3 Binarity 

One intriguing possibility is that we are seeing variability in- 
duced by a companion orbiting Kelu-1. A companion could 
modify Kelu-ls brightness in two ways; either by partial or 
total eclipsing, or via the gravitational pertubation of Kelu- 
I's photosphere, causing ellipsoidal variability (see Hilditch 
2001). We have not been able to derive any orbital con- 
figuration which could reproduce the observed lightcurve 
amplitude and shape, and we therefore reject an eclipsing 
companion. In the case of ellipsodial variability, we would 
observed two maxima per rotation, and hence the orbital 
period would be 3.6 hours. It is reasonable to assume the 
orbits will be circular due to tidal effects, and we can there- 
fore directly calculate the orbital seperation as; 



a = 0.48J? 



■0 



Ml 



O.O65M0 



1/3 



(s.eh) 
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where M\ is the primary (Kelu-1) mass, and P is the orbital 
period. 

Unfortunately, when fitting for ellipsodial variability, 
there is a degeneracy between the mass of the secondary 
object and the orbital inclination of the system. If, however, 
we assume the secondary companion is also a brown dwarf 
(lack of X-ray flux rejects the possibility of a massive dark 
companion such as a neutron star or black-hole), we can 
place limits on its mass. Figure ^ shows the possible sec- 
ondary mass as a function of primary mass. Companions are 
excluded from the hatched region for two different reasons: 
(1) if M\ > 45A4jyp, the minimum mass is determined by 
the need to produce 1.1% ellipsodial variability without the 
secondary causing eclipses (i<65°), (2) for Mi < 45Afjyp, 
brown dwarfs cannot exist in the hatched region, as they 
would overflow their Roche-lobes, transfering matter to the 
primary. It is not clear what the final outcome of such an 
interaction would be, but the lack of activity indicates this 
cannot currently be happening. An object more compact 
than a brown dwarf or gas-giant planet could survive in 
this region. We note that these limits are based on sim- 
plistic models and are very sensitive to the exact amplitude 
of Kelu-1 's variability. They do however indicate that any 
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Figure 9. Allowed masses for a companion to Kelu-1. If Kelu-1 
has a mass greater than ~45Mjup, the minimum mass limit is set 
by the need to induce 1.1% elliposidal variations with i<65° . How- 
ever, if Kelu-1 has a mass less than ~45Mj„p, brown dwarfs/gas- 
giant planets arc excluded from the hatched region as they would 
overflow their Roche-lobe. 
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companion to Kclu-1 would be a brown dwarf rather than 
massive planet (i.e. M2 > ISMj^p). 

Pure elUpsodiaJ variabihty produces a sinusoidal 
hghtcurvo, which is statistically consistent with Kelu-l's 
folded lightcurve. If a secondary companion is responsible 
for the variability, future observations will return the same 
amplitude and ephemeris. We note that a binary system 
would have an orbital velocity of ~130kms~^, and would be 
capable of mimicing the observed 60kms~^ line broadening 
over a 1 hour integration. 

4.4 Surface structures 

Doppler imaging codes have been used extremely successful 
to mapping and understand the surface features on solar- 
type low-mass stars. However, the process requires very high 
quality photometry and spectroscopy, combined with a full 
understanding of photospheric physics; all of which seem a 
long way off for brown dwarfs. 

Although mapping the surface of a brown dwarf is a 
distant goal, we can make some headway. By synthesising 
photometry from theoretical spectra (Allard et al. 2001) of 
two limiting cases; 1) dust is suspended in the atmosphere 
(DUSTY), and 2) dust forms, but immediately settles below 
the photosphere (COND), we can investigate what degree 
of surface structure is required to reproduce the observed 
variability. 

Kelu-l's spectrum is best matched by a iJf/iS'T'y atmo- 
spheric model (Baraffe et al. 1998), so we assume a model 
in which variability is caused by clear ( COND) "holes" in 
a photosphere dominated by DUSTY "clouds". The COND 
spectrum is 0.55 mag brighter than the DC/^Ty spectrum 
through our filter. The amplitude of 0.012 mag we observe in 
Kelu-1 corresponds to COND "holes" covering ~1.7% of an 
otherwise I>[/STy photosphere. This value is a lower limit, 
as less perfect clear patches would have a smaller effect on 
the variability, requiring them to cover larger fractions of 
the surface. 



Note that the covering fractions above actually repre- 
sent the difference in covering fraction between maximum 
and minimum light. They are therefore only representative 
of surface features on ~180° scale. It may be that the sur- 
face has a far more significant small scale structure, which 
simply averages out over large regions. 

4.5 The radius of Kelu-1 

Interpreting the l.Sh variability as the rotation period al- 
lows us to make one of the first firm tests of brown dwarf 
evolutionary theory. Evolutionary models (e.g. Chabrier et 
al. 2000) predict that after ~100Myr, brown dwarfs (inde- 
pendant of mass) reach a stable radius of ~O.1R0. We can 
combine the rotation period and the rotational velocity to 
produce an observational lower limit (due to the unknown 
inclination) on Kelu-l's radius. For the observed values of 
1.80±0.04h and 60±5kms"^ (Basri et al. 2000), Kelu-1 must 
have R > 0.09 ±0.0171©. The evolutionary models are there- 
fore in agreement with our observations of Kelu-1. This is 
not the case for the observations of BRI0021-0214 made by 
Martm et al. (2001). They require a radius of at least O.MR© 
to fit their observations with meta^stable photospheric sur- 
face features. We note that this argument would require the 
features responsible for photometric varibility to have the 
same rotation period as the photosphere. 

Alternatively, we can invert this argument and, assum- 
ing a model radius of O.1±O.O5R0 (to include a range of pos- 
sible ages and masses), determine Kelu-l's inclination to the 
line of sight as 53 < i < 77. This means that between 80% 
and 97% of the surface can cause variability. The remain- 
ing fraction (polar regions) is either permanently in view, or 
never visible. 



5 CONCLUSIONS 

It is now clear that variability from brown dwarfs can be 
detected, but that photometry of better than 1% is required 
to do so. In the near future, futher studies of variability 
and rotation in brown dwarfs should greatly increase our 
understanding of their physics. 

We have detected a strong periodicity of 1.80±0.05 
hours in differential photometry of the L2 brown dwarf Kelu- 
1. We have investigated four possible mechanisms to explain 
this variability: 

(i) surface inhomogeneity moderated by meteorology and 
variable dust formation; 

(ii) surface inhomogeneity moderated by magnetic 

starspots; 

(iii) light curve variability due to gravitational distortion 

of Kelu-l's envelope by a close companion; and 

(iv) light curve variability due to an eclipsing binary 

Mechanisms (ii) & (iv) seem unlikely explanations, but 
we are unable to concusively differentiate between mech- 
anisms (i) & (iii). Ellipsodial variability, mechanism (iii), 
would produce a twin peaked lightcurve, giving a period of 
3.6±0.1 hours. This mechanism will give stable and repeat- 
able photometric variability in future epochs. Alternatively, 
mechanism (i) associates the 1.80±0.05 hour period with 
the rotation period, which is consistent with the rotational 
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velocity of 60 kms^^ and theoretical radius of ~0.1Rq, in- 
dicating an inclination in the range 53° < i < 77° . 

Over the duration of our observations, the general shape 
and period of the lightcurve are unchanged (at least to 
within the measurement noise), implying the process caus- 
ing the modulations is also stable on this timescale. This 
is in contrast to previous variability observations where no 
periodicity, or a period that changes on the order of the 
observation length, has been measured. 

The two explanations wo have presented lead to differ- 
ent predictions for future observations of Kelu-1. Variabil- 
ity induced by a secondary companion will be completely 
repeatable at future epochs, whereas long term evolution 
of surface features will result in a secular changes in the 
lightcurve. Further observations of Kelu-1 will be a powerful 
discriminant between the two hypotheses we have presented. 
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